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Abstract 

We calculate the 0{as) QCD and 0{acw) electroweak one-loop corrections in the Standard 
Model framework, to the production of an intermediate Higgs boson associated with tt pair via 
77 fusion at an electron-positron linear collider (LC). We find the 0{as) QCD corrections can 
be larger than the 0{aew) electroweak ones, with the variations of the Higgs boson mass mh 
and e'^e" colliding energy ^/s. Both corrections may significantly decrease or increase the Born 
cross section. The numerical results show that the relative corrections from QCD to the process 
e+e^ ^ 77 — > tth'^ may reach 34.8%, when ^/s ~ 800 GeV and nih = 200 GeV, while those from 
electroweak can be -13.1%, -15.8% and -12.0%, at = 800 GeV, 1 TeV and 2 TeV respectively. 
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I Introduction 

After the discovery of top quark, directly searching for Higgs boson and studying its property are 
the main goals of the high energy colliders. In the electroweak minimal standard model(MSM), 
Higgs mechanism generates electroweak symmetry breaking and the Yukawa coupling terms be- 
tween Higgs boson and fermions in the Lagrangian H]- In the SM, the Yukawa term reads 
Cy = — Ylf '^Z (1 + -^(^)/^) ''Pf'^f- There the coupling strength of the fermion-Higgs Yukawa cou- 
pling f — f — is predicted as g^jf^ = rrif/v at the tree level, where v = {y/2GF)^^^'^ — 246 GeV is 
the vacuum expectation value of the Higgs field. Since the top quark is the heaviest fermion, the top 
quark Yukawa coupling gr^f/j should be the strongest one among all the fermion-Higgs couplings, e.g., 
9tth ~ compared for example with g^j^^ ~ 4 x 10^^. Therefore, the Higgs boson production 

associated with a top-quark pair production process is particularly important in collider physics for 
probing the top-quark Yukawa coupling with the intermediate mass of Higgs boson. 

Recently, LEP2 experiments have provided a lower bound of 114.4 GeV for the SM Higgs boson 
mass at the 95% confidence leveljS]. The future linear colliders (LC) will continue the work in searching 
for Higgs boson and studying its property. There have been already some detailed designs of linear 
colliders, such as NLCjl], JLC|n|, TESLAjHl and CLIC[3. Although the cross section for e+e~ tthP 
process is small at a LC, about 1 fb for ^/s = 500 GeV and = 100 GeV jHHHEl) it has a distinctive 
experimental signature and can potentially be used to measure the top-quark Yukawa coupling with 
intermediate Higgs mass at a LC with very high luminosity. 

As we know the apparently the clean signal for light Higgs boson production associated with a 
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top quark pair is e^e~ tih? bbbbW^W~ in both semileptonic and fully hadronic decay channels. 
This leads to multi-jet event topologies involving at lest 6 or more jets in the final state, with > 4 
b-jets and multi-jet invariant mass constraints. Obviously, its measurement has many difficulties, 
for example, the tiny signal with backgrounds about 3 orders of magnitude larger, the limitations of 
jet-clustering algorithm in properly reconstructing multi-jets in the final state, and the degradation 
of b-tagging performance due to hard gluon radiation and jet mixing. The potential backgrounds are 
from the ttZ^ and tt productions. The dominant electroweak background is: 

77 ^ ttZ^ bbW+W~Z° bbbbl^uqq'. (1.1) 

The largest background is from radiative top quark decay: 

-f-f^tt^ bbW^W^g bbbbl^uqq'. (1.2) 

Since the b-jets resulting from the gluon splitting are logarithmically enhanced at low energy, cuts 
on the jet energy are efficient at eliminating this background The weak corrections to 77 tihP 
should involve the real emissions of W and gauge bosons. These may be other sources of electroweak 
backgrounds. Our calculation shows that the cross section of the real emission of Z^ gauge boson in 
tih? production via 77 collision is under 0.1% of the 77 tih^ process, and can be neglected in our 
analysis, while for the real W emission tih^ production the cross section is in the same order as that of 
the 77 — > tth^ process. In principle, these backgrounds can easily be removed by using the constraints 
due to the W, Z'^,t and /i° masses with perfect b-tagging and reconstruction of multi-jets. However, 
this may not be true in practice. The experimental situation of the signal detection of tih^ production 
will be even worse, considering the following two fields. Firstly, the detectors have a finite coverage in 
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polar angle. Some of these real W bosons can be emitted so close to the beam direction that part of 
their decay products may well escape detection. Secondly, all three particles in the tth^ final state are 
highly unstable, and would decay in the detector. If rrih ^ 130 GeV, it would mostly produce bb pairs. 
If above, it would yield four fermions, either leptons/neutrinos or quarks, principally via — > WW* 
decays. Besides, top quark(anti-quark) decays into three objects, a 6-quark and the decay products 
of the W, leptons/neutrinos or quarks. Hence, at tree level in these decays, we are looking at 8- or 
10-fermion final states. In the most likely case in which the fermions arising from the W decays are 
quarks, additional (mainly gluon) radiation will take place, so that at detector level more than 10 
jets may be extracted. Therefore, the additional real W radiation and ttZ^ production entering the 
detector region could be uneasy to resolved through their decay products in the data samples of tth^ 
events. 

Expected experimental accuracy for determination of the tih^ coupling in e~^e~ — > tth^ process 
has been discussed in many literatures for specialized linear collider. They show that the experimental 
accuracy depends very much on the b-tagging efficiency. For example, the top quark Yukawa coupling 
can be measured to 6 — 8% accuracy with integral luminosity 1000 fb^^ at an e~^e~ linear collider with 
^/s = 1 TeV, assuming 100% efficiency for b-jet tagging and including statistic but not systematic 
errors. The accuracy of the measurement drops to 17 — 22% if only a 60% efficiency for b-tagging is 
achieved jlll I12j . The references |131 1141 [T3] also stated that the precise determination of the top- 
Higgs coupling via the measurement of e'^e~ — > tih^ process can reach the accuracy of few percent. 
Therefore, if we assume the machine and detector are very efficient and the background processes 
can be distinguished substantially, the evaluation of radiative corrections could be significant for the 
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accurate experimental measurements of e^e~ — > tth^ process. In Ref-I^], S. Dawson and L. Reina 
presented the NLO QCD corrections to process e~^e~ — > tt/i" . And in references |161 1171 118j the 
electroweak corrections to the process e^e~ — > tih^ are calculated. The supersymmetric electroweak 
corrections to process e^e~ — > til^ were already discussed by X.H. Wu, et al[19j. 

An e^e~ LC can also be designed to operate as a 77 collider. This is achieved by using Compton 
backscattered photons in the scattering of intense laser photons on the initial e^e~ beams. The 
resulting 7 — 7 center of mass system (CMS) energy is peaked at about 0.8-^5 for the appropriate 
choices of machine parameters. In 77 collision mode at the high energy peak, we may get approximately 
the same luminosity as that of e^e~ collision. With the new possibility of 77 collisions at electron- 
positron linear colliders, the production process e^e~ ^ 77 ^ tt/i'^ offers another approach to probe 
directly the top-Higgs coupling in addition to e^e~ — > tt/i^ andpp(pp) — > ttlrS^ processes. To detect the 
Higgs boson associated with a top-quark pair in high colliding energy, 77 collision has an outstanding 
advantage over e^e' collision due to its relative larger production rate. The reason is that at the 
tree level the e^e~ — > ttl^ process has a 's-channel suppression' from the virtual photon and 
propagators, especially for the heavy masses of the final particles. Therefore, we can conclude that 
the e^e~ ^ 77 ^ ttl^ process provides a better approach than e^e^ collision to produce tt/i^. 
Similar with the measurement of the process e^e~ — > tt/i'^ , the evaluation of radiative corrections for 
process 77 tth^ is also significant for the accurate experimental measurements of top quark Yukawa 
coupling. 

The Born cross section of e~^e~ ^ 77 — > tth^ process was calculated already in previous work of 
Ref ■ |2Uj ■ In this paper we neglect the real gauge boson radiation effects and present the calculations 
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of the 0{as) QCD and 0{aew) electroweak corrections to e~^e~ ^ 77 ^ tth^ in the SM. We draw a 
comparison between our numerical results of Born cross sections of 77 — > tth^ and those in Ref . j2Uj . 
The paper is organized as follow. In section 2, we present the notations and analytical calculations 
of the QCD and the electroweak radiative corrections. The numerical results and discussions are 
presented in section 3. Our conclusions are given in section 4. The numerical comparison of the Born 
cross section of process e+e" — > 77 ^ tth^ are presented in Appendix. 

II Analytical Calculation 

II. 1 Calculations of the lowest order of the 77 tth? subprocess 

The subprocess 77 — > tihP at the lowest level occurs through the u- and t-channel mechanisms involving 
Higgs boson bremsstrahlungs originated from different positions on top-quark lines. The tree level 
diagrams are drawn in Fig^ but the corresponding diagrams with interchange of the two incoming 
photons are not shown. 




Figure 1: The lowest order diagrams for the 77 — > tth^ subprocess. 
We denote the subprocess 77 tthP as 

liPl) + l{P2) ^ t{ki) + i{k2) + h\k^). (2.1) 

The four-momenta of incoming electron and positron are denoted as pi and respectively, and the 
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four-momenta of outgoing top-quark, anti-top-quark and Higgs boson are represented as fci, k2 and 
fes correspondingly. All these momenta obey the on-shell equations Pi = P2 = 0, /cf = /c^ = 

kl = ml 

The amplitudes of the corresponding t-channel Feynman diagrams (shown in Fig^a-c)) of the 
subprocess 77 — > tih^ are represented by 

j^ia) ^ _ e^Qlmt 1 1 



2mw sin 6w {ki — piY — rrq {p2 — ^2)^ — 
xu{ki)i{pi){^i + mt){^2 -h + mt)i{p2)v{k2), (2.2) 



j^{b) ^ e^Qlmt 



2mw sin^vK (^1 — PiY — (^2 + k^Y — ml 
xu{ki)i{pi){]/^i + mt)i{p2){-h -h + mt)v{k2), (2.3) 



2mw svnOw {ki + k^Y - ml {p2 - ^2)^ - 
xu{ki){h + h + mt)i{pi){h -h + mtW2)v{k2), (2.4) 

where Qt = 2/3 and the corresponding amplitudes of the u-channel Feynman diagrams of the subpro- 
cess 77 — > tth^ can be obtained by the following interchanges. 

M'^^ = M^^\pi ^ P2), A^l") = Mf\pi ^ P2), M^:'> = M^ipi ^ P2). (2.5) 

The total amplitude at the lowest order is the summation of the above amplitudes. 

c t 



i=a,b ]=u 

Although the previous calculations for the 77 — > tih? subprocess and e+e~ ^ 77 — > tih? process at 



the lowest order were presented by Kingman Cheung[2n], we made the numerical comparison with his 



results yet. We calculated the Born cross section of the process 77 tihP by using Feynman gauge 
and unitary gauge to check the gauge invariance, and adopting FeynArts 3,21^ and CompHEP|22j 
packages respectively. We found our results are in good agreement with each other, but not coincident 
with Cheung's. The numerical comparisons are presented in Appendix Table 3. 



II. 2 Calculations of the 0{as) QCD corrections of the 77 tth!^ subprocess 




Figure 2: The QCD pentagon diagrams for the 77 tthP subprocess, whose amplitudes include five- 
point tensor integrals of rank 4. The corresponding diagrams with interchange of the two incoming 
photons are not shown. 

The 0{as) QCD one- loop Feynman diagrams of the subprocess 77 tthP are generated by 
FeynArts 3 |21j . There are 84 Feynman diagrams with 0{as) corrections of the virtual one-loop QCD 
corrections, which involves the vertex correction, internal propagator self-energy correction, box and 
pentagon diagrams. The Feynman graphs which generate amplitudes including five-point integrals of 
rank 4 are shown in Fig|2 as a representative selection. The amplitude of the subprocess 77 — > tthP 
including virtual QCD corrections to 0{as) can be expressed as 



(2.7) 



Mqcd = Mo + -^CfM^c^cd- 
where Cp = 4/3, The term j^CfMqqj^ is the amplitude contributed by the QCD one-loop Feynman 



diagrams and the QCD renormalizations of top-quark wave function, mass and t — t — hP Yukawa 
coupling. We define the relevant renormalization constants as 



m,,o = mt + (5m,(,), ^ = (1 + = (1 + = 1^ + ^-11^ . 



2 t(9)) 



V V 



(2.8) 



With the on-mall-shell renormalized condition we get the 0(as) QCD contributed parts of the renor- 
malization constants as 



Sm 



t{9) 



mt 



Re (Sf(^)(mf ) + S,^^)(m?) + 2Sf(^^ 



(mf )) , 



(2.9) 



6Z, 



t{9) 



d 



d 



\p'^=m1 ■> 



\p'^=m^i 



(2.10) 



(2.11) 



Re takes the real part of the loop integrals appearing in the self energies only. Here we define the 
renormalized top-quark irreducible two-point function as 



(2.12) 



where a and (3 are the color indices of the top quarks on the two sides of the self-energy diagram, 
Pl,r = (1 T 75)/2- The unrenormalized top-quark self energy parts contributed by 0{as) QCD read 



Sf(^)(p2) = S,^,)(p') = (-1 + 2i?o[p,0,mt] + 25i[p,0,mi]) 



(2.13) 



and 



|^(l-2Sob,0,mt]) 



(2.14) 



The corresponding contribution part to the cross section at 0{as) order can be written as 

= ^oe" = I d^sEl^e {MlM^hn) , (2.15) 

I "1 1 spin 

where d$3 is the three-body phase space element. The bar over summation recaUs averaging over 
initial spins. 

The virtual QCD corrections contain both ultraviolet (UV) and infrared (IR) divergences in general. 
To regularize the UV divergences in loop integrals, we adopt the dimensional regularization in which 
the dimensions of spinor and spacetime manifolds are extended to D = 4 — 2e. We have verified the 
cancellation of the UV both analytically and numerically. Then we get a UV finite amplitude including 
0{as) virtual radiative corrections. 

The IR divergence in the ■M.qqjj of the process 77 tihP is originated from virtual gluon cor- 
rections. It can be exactly cancelled by including the real gluon bremsstrahlung corrections to this 
subprocess in the soft gluon limit. The real gluon emission process is denoted as 

l{pi) + l{V2) ^ t{ki) + t{k2) + h^h) + g{k), (2.16) 

where the real gluon radiates from the internal or external top(anti-top) quark line, and can be 
classified into two parts which behave soft and hard natures respectively. In order to isolate the soft 
gluon emission singularity in the real gluon emission process, we adopt the general phase-space-slicing 
method j33j . in which the bremsstrahlung phase space is divided into singular and non-singular regions. 
The cross section of the real gluon emission process (|2.16() is decomposed into soft and hard terms 

= AaZf + A^£;- = MsZf + C-). (2.17) 
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We adopt the soft gluon approximation method with a cut IS.E in numerical calculations for the soft 
emission corrections (A;" < IS.E). With this approach the real gluon radiation from internal color lines 
does not lead to IR singularities and can be neglected in this approach. We find the contribution of 
the soft gluon emission process is [SH 



ki 



ko 



-\ 2 



ki ■ k k2 ■ k 



(2.18) 



in which AE is the energy cutoff of the soft gluon and k^ < AE <^ y/s. k^ = \J\k\'^ + rrig is the 
gluon energy. Here we have introduced a small gluon mass nig to regulate the infrared divergences 
occurring in the soft emission. The integral over the soft gluon phase space have been implemented, 
the analytical result of the soft gluon corrections to 77 — > tih^g is presented as [5] 

xki ■ /C2 



-2 log 



27r 



'm„ 



1 



m^{x'^ — 1) 



2 1 I ^^2 I '^2 1 



M V^2 + l^2 



+ 



+ Li2 1 



u — \u\ 



Axki ■ k2 
m^(x^ — 1) 

u=xk\ '\ 

\ ' 

u=xk2 ) 



u — \u\ 



vP + \u\ 





fe? 




1^1 1 


+ Li2 fi- 



u + n 



where 



and X is the solution of 



mf (x^ — 1) 

2(x/fc?-A;0)' 



mfix'^ + 1) - 2xki ■k2 = 0, 
which should satisfy the constraint of 



"-2 



> 0. 



(2.19) 



(2.20) 



(2.21) 
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We checked numerically the cancellation of IR divergencies and verified that the contribution of 
these soft gluon bremsstrahlung corrections leads to an IR finite cross section which is independent 
of the infinitesimal gluon mass mg. The hard gluon emission cross section Acj^^J^ for Eg > AE, is 
calculated numerically by using Monte Carlo method. The statistic error is controlled under 0.3%. 

Finally the UV and IR finite total cross section of the subprocess 77 — > tth^ including the 0{as) 
QCD corrections reads 

= ao + Aa^^^ = ao + Aa^^^ + Aa.^^^ = ao(l + (2.22) 
where J'^'^^ = 6^^f^ + (5^^^ + ^^^^ is the QCD relative correction of order 0{as)- 

II. 3 The calculation of the C(q!c„) one-loop corrections to the 77 — > tth'^ subprocess 

In addition to the QCD corrections, we also calculate the 0{aew) one-loop electroweak corrections to 
the subprocess 77 tihP . We use again the package FeynArts 3j2J| to generate the electroweak one- 
loop Feynman diagrams and the relevant amplitudes of the subprocess 77 — > tth^ automatically. The 
electroweak one-loop Fenyman diagrams can be classified into self-energy, triangle, box and pentagon 
diagrams. The pentagon diagrams, whose corresponding amplitudes include five-point tensor integrals 
of rank 4, are depicted in Fig|31as a representative selection. In our electroweak correction calculation 
we use also the t'Hooft-Feynman gauge and adopt the definitions of one- loop integral functions of 
Ref . j29j ■ The one-loop level virtual electroweak corrections to 77 — > tth^ can be expressed as 

A^.T = ao6^^ = [ dcI>3^Re (mIm^-^) , (2.23) 
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where pi is the c.m.s. three-momentum of one of the incoming photons, d$3 is the three-body 
phase space element, and the bar over summation recalls averaging over initial spins. -M.^^ is the 
amplitude of the electroweak one-loop Feynman diagrams, including self-energy, vertex, box, pentagon 
and counterterm diagrams. 

Analogously to the case of the QCD correction calculations shown in last subsection, we adopt 
the dimensional regularization scheme to regularize the UV divergences in loop integrals. We assume 
that there is no quark mixing, the KM-matrix is identity matrix and use the complete on-mass-shell 
(COMS) renormalization scheme |32]) in which the electric charge of electron e and the physical 
masses m^y, mz, ruh, rrit et al., are chosen to be the renormalized parameters. The relevant field 
renormalization constants are defined as 



222 

eo = {l + 5Ze)e, mf^Q = mh + Smf^, mtfl = nit + 5mt, myirQ = 'miY + 6mw, 
fn'zfl = "^1 + ^0 = ]^5ZazZ + (1 + ^6Zaa)A, /iq = (1 + ^6Zh)h, 

= (1 + \5Z^)t\ t§ = (1 + ^6Z^)t^. (2.24) 

With the on- mass-shell conditions, we can obtain the renormalized constants expressed as 

5mlr = ReJ:^iml,), 5ml = ReT,^^ {ml) , 5Zaa = - ^^^^^^""^ 



gp2 IP^=0' 

6Zzz = -Re ^^ff\ ,.=^., 6ZzA = 2^i^, 5 Z az = -2Re^^^ , 



z 

5mh = Re^\ml), 5Zh = -Re—^^\p2=^2. (2.25) 



The renormalization constants of the wave function and mass of top-quark can be evaluated from 
Eas. (|2.9|) - (|2.11j) upon replacing the QCD top-quark self-energies(Il^^^-|, S^^^ and 5]^^^^ by the elec- 

13 



troweak corresponding ones (Sj', and Sf ), respectively. And the explicit expressions of the elec- 
troweak self energies in the SM concerned in our calculations can be found in the Appendix B of 
Ref ■ j32j ■ The UV divergence appearing from the one-loop diagrams is cancelled by the contributions 
of the counterterm diagrams in our calculations. We have verified both analytically and numerically 
that the final cross sections including 0{aew) virtual radiative corrections and the corresponding 
counterterm contributions are UV finite. The further verification of the correctness for the one-loop 
calculation is made by probing the gauge independence of cross section via changing the value of ^ in 
gauge, and the results are coincident with each other very well. 

Analogous to the calculation of the QCD corrections, the IR divergence in the subprocess 77 
tth^ originates from virtual photonic corrections is cancelled by the real photonic bremsstrahlung 
corrections in the soft photon limit. We use also the general phase-space-slicing method |33j and 
divide the phase space into singular and non-singular regions. Then the cross section of the real 
photon emission subprocess (77 — > tth^'y) is decomposed into soft and hard parts. 

Ao-real = A(7soft + A<Thard = ^oi^]oh + ^hard)" (2.26) 

By using the soft photon(A;'^ < AE) approximation, we get the contribution of the soft photon emission 
subprocess expressed as jSH EU 

ki k2 



soft - -"C'O „ 2 / ^ ^ 

J\k\<AE 2ko 



ki- k k2 ■ k 



(2.27) 



in which AE is the energy cutoff of the soft photon and k^ < AE ^ Qt = 2/3 is the electric 



charges of top quark, k^ = ylkl"^ + rrii^ is the photon energy. Therefore, after the integration over the 
soft photon phase space, we obtain the analytical result of the soft corrections to 77 tthP . Actually, 

14 



the real photonic emission correction can be deduced from the real gluon emission corrections upon 
replacing the factor CpOs by Q'^Oew The cancellation of IR divergencies is verified and the results of 
the cross section show the independence on the infinitesimal photon mass in our calculation. 

Since sometimes the QED radiative contributions can be quite large, the investigation of the 
genuine weak corrections quantitatively would help us to understand the origination of the large 
electroweak corrections. The QED correction to the subprocess 77 — > tth^ is gauge invariant and 
comprises three parts: (1) photonic virtual radiations of final and internal top-quarks, (2) real photonic 
radiations from final and internal top-quarks, (3) the interference of final and internal real photon 
radiations. Actually, the total ©(aew) order QED corrections for the subprocess 77 tih^ , can 
be obtained numerically from the results of the 0{as) order QCD corrections as presented in section 
II. 2 through multiplying a factor of ^"(J^y^^ ■ Then the genuine weak corrections can be evaluated 
by subtracting the QED corrections from the 0{aew) QED corrections. We define the genuine weak 
relative correction as, 

L = S- 6^^^ = 6 - - - Crl (2.28) 

where S^^f^ is the relative correction contributed by the QED one-loop diagrams including virtual 
photon exchange and the corresponding parts of the counter terms. (J^,//^ ^i^d S^^^j^ are the relative 
corrections of the soft and hard real photon emissions, respectively. 

II. 4 Calculations of process e~^e~ — > 77 ^ tih^ 

By using the laser back-scattering technique on electron beam, an e'^e'^ LC which has c.m.s. energy 
of hundred GeV to several TeV can be transformed to be a photon collider |241 1251 [26] . By integrating 
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over the photon luminosity in an e~^e hnear cohider, the total cross section of the process e'^e — > 
77 tth^ can be obtained in the form as 

a{s) = / dz— -^(t(77 tth at s = z'^s) (2.29) 

where £^0 = +m/i, and ^/s{\/^) being the e'^e~ (77) c.m.s. energy. ^^^^^ is the distribution function 
of photon luminosity, which is defined as: 

^ = 2z ^F,/,(x)F,/,(zV^) (2.30) 

For the initial unpolarized electrons and laser photon beams, the energy spectrum of the back scattered 
photon is given by |77] 

1 



D{0 



1 Ax Ax^ 
1 -x+- -7- r + 



(2.31) 



1-X i{l-x) C2(l-X)2_ 

where x = 2uj/y/s is the fraction of the energy of the incident electron carried by the back-scattered 
photon, the maximum fraction of energy carried by the back-scattered photon is Xmax = 
e/(l+e), and 



m) = (1 - 7 - 7^)in(i +e) + + 7 - .... . .,2 ^ (2-32) 



e = (2.33) 

rrie and \/i/2 are the mass and energy of the electron, ujq is the laser-photon energy. In our evaluation, 
we choose loq such that it maximizes the backscattered photon energy without spoiling the luminosity 
through e~^e~ pair creation. Then we have ^ = 2(1 + V2), 

Xmax — 0.83, and D{(,) ~ 1.84, as used in 

Ref.HH]- 
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Ill Numerical results and discussions 



For the numerical calculation we use the following input parameters |35| 

aew(O)"^ = 137.03599976, mw = 80.423 GeV, mz = 91.1876 GeV, 

me = 0.510998902 MeV, = 105.658357 MeV, = 1.77699 GeV, 

m„ = 66 MeV, = 1.2 GeV, rrit = 174.3 GeV, 

rud = 66 MeV, = 150 MeV, rrib = 4.3 GeV, 



as[m 



I) = 0.117186. 



where we use the effective values of the light quark masses (m„ and m^) which can reproduce the 
hadron contribution to the shift in the fine structure constant aew("^^) [SHI) and take the pole mass 
of top quark(m( = 174.3 GeV) to determine the tth^ coupling strength. The QCD renormalization 
scale fi is taken to be {2mt + mh)/2 and the running of the strong coupling as(/i^) is evaluated at the 
two-loop level(M S scheme) with five active flavors. 

The numerical results of the cross sections with QCD and one-loop electroweak radiative corrections 
for the subprocess 77 — > tih^ are plotted in Fig^and Fig|SJ respectively. The full, dashed, and dash- 
dotted curves correspond to the cases with rrih = 115, 150 and 200 GeV, correspondingly, and 77 
colliding energy \/I runs from the value little larger than the threshold (2mj + rrih) to 1.8 TeV. For 
each line type there are two curves, the upper curve (in the region \/I > 1 TeV) is for the Born 
cross section and the lower one for the QCD corrected cross section. As indicated in Fig01 the QCD 
corrections can increase (when \f§ < 650 GeV) or decrease the cross sections of subprocess 77 tih^ 
(when \f§ > 900 GeV), while FigO shows that the one-loop electroweak radiative corrections always 
reduce the Born cross sections in the plotted energy range of The curves in both FigHand FigEl 
show that all the Born, QCD and electroweak corrected cross sections decrease with the increment of 
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the mass of Higgs boson mh. The curves for mh = 115 GeV increase rapidly to their corresponding 
maximal cross section values, when the 77 colliding energy \/I goes from the threshold value to the 
corresponding position of peak. The curves for ruh = 150 GeV have platforms when \/I are larger 
than 1000 GeV, and for = 200 GeV both Born and one-loop corrected cross sections increase 
slowly in the whole plotted range of 

We define the QCD relative corrections to the subprocess 77 — > tth^ and parent process e^e~ — > 
77 tth^ as 

-QCD i. ^QCD 
^QCD ^ ^QCD ^ -^0 ^ (3.1) 

respectively, and the electroweak relative corrections to the subprocess 77 tihP and process e^e" — > 
77 tihP as 

J (T - ao cj - fJo 

= — , = , 

Co cfq 

(3.2) 

correspongdingly. The 0{as) QCD relative corrections and ©(oew) electroweak relative corrections 
to the cross sections for 77 tih? subprocess, corresponding to Fig^ and FigO are depicted in 
Fig El and Fig|7fa-b), respectively. We can read from FigEl that the QCD relative corrections to 
the subprocess 77 — > tih^ decrease from 111% to —13.8% when the c.m.s. energy \/I increases from 
the threshold energy to 1.8 TeV. From Fig. 7(a) we can see that in the plotted colliding energy range, 
the O(aew) order electroweak relative corrections to subprocess 77 tih? can reach —7.98% and 
— 16.5% for rrih = 115 GeV and 200 GeV, respectively. The maximal electroweak absolute relative 
corrections to the cross sections |(5|max and the corresponding ^fI positions for subprocess 77 tih^ 
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with ruh = 115, 130, 150, 170, 200 GeV are listed in Tabled In Fig.7(b) the ©(oew) order QED and 
genuine weak relative corrections are depicted, respectively. It shows that the 0{aew) QED relative 
corrections are very small comparing with the genuine weak relative corrections, and can only reach 
1.28% for the curve of nih = 200 GeV at the position of \/I ~ 580 GeV. 

We also depicted the QCD and electroweak relative corrections to the cross sections of subprocess 
77 tih? as the functions of the mass of Higgs boson rrih with \/I = 500, 800, 1000, 2000 GeV in FiglHl 
and FigEl respectively. Both curves for \/I = 500 GeV in FiglHl and FiglHl are truncated before the 
position of ~ 150 GeV, because the channel 77 tth^ cannot be opened when \/I < irifi + 2mt. 
On the curves for \/I = 800, 1000, 2000 GeV in FiglHl there exist two resonance peaks on each 
curve at the positions nih ~ 2m z and m/j ~ 2m\Y. Since we didn't consider the widths of and 

gauge bosons in loop calculation, the genuine one-loop weak corrections in the vicinities of the 
thresholds at nih = 2myY and m/j = 2mz shown in FiglHl are untrustworthy. In FigEl the QCD 
relative corrections for \fl = 500 GeV are rather large, and can reach the value larger than 110% 
when the ruh = 140 GeV and \/l = 500 GeV^, while the QCD relative corrections are in the range 
1.74% ~ -14.0%, when \fl = 800, 1000, 2000 GeV. FigH shows that the curves of the electroweak 
relative corrections for \/I = 800 GeV and 1000 GeV go down from —6.48% to —16.0% and from 
—7.01% to —13.9% respectively, when rrih varies from 100 GeV to 200 GeV. While the electroweak 
relative corrections for ^/I = 2000 GeV are relative stable( about minus a few percent) except in the 
Higgs mass regions satisfying resonance conditions. 

Fig llOl and Fig llll show the cross sections including 0{as) QCD and 0{aew) radiative corrections 
and for e'^e~ ^ 77 ^ tih^ process versus the e^e~ colliding energy y^, respectively. Both figures 
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demonstrate that the Born and radiative corrected cross sections increase with the increment of ^/s, 
and the 0{as) QCD radiative corrections for different values of nih can reduce or increase the Born 
cross sections, but the 0{aew) electroweak corrections reduce the Born cross section only in the plotted 
range of electron-positron c.m.s. energy ^/s. The QCD and electroweak relative corrections for the 
e~'~e~ — > 77 ^ tth^ parent process corresponding to Fig. 10 and Fig. 11, are depicted in Fig. 12 and 
Fig. 13, respectively. We can see that the QCD relative correction can be very large near the threshold 
colliding energy, while the electroweak relative corrections can reach —15.9% for nih = 200 GeV in the 
vicinity of y/s ~ 1 TeV. The maximum absolute electroweak relative corrections for different Higgs 
boson mass values in the plotted colliding c.m.s. energy range can be read out from this figure and 
are listed in Table [21 

The QCD and electroweak relative corrections to the cross section of the process e~^e~ ^ 77 — > tih^ 
as the functions of the Higgs boson mass rrih are depicted in Fig il4l and Fig ll5l with y/s = 800, 
1000, 2000 GeV, respectively. From Fig |14l we can see that the QCD one-loop relative correction for 
^/s = 800 GeV can reach 34.8% at the position of m/j = 200 GeV. In Figil5ltwo resonance peaks appear 
again on each curve at the positions of rrih ~ 2mz and m/j ~ 2mw because of the resonance effects, 
and the correction values in these vicinities are untrustworthy. For the curves of nih = 800 GeV and 
nih = 1000 GeV, the relative correction decrease from —2.66% and —4.51% to —13.1% and —15.9%, 
respectively, when nih varies from 100 GeV to 200 GeV. The relative correction for ^/s = 2000 GeV 
varies in the range of [—6.15%, —12.0%] when rrih goes from 100 GeV to 200 GeV, except in the 
mass regions satisfying the resonance effect conditions. 

Scale dependence of the i^-factors (= "'^^ ) of the total cross sections for the process e~^e~ — > 
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77 tth^ at LC {nih = 115 GeV) is plotted in Fig |16l The full-line, dashed-line and dash-dotted-line 
correspond to ^/s = 800 GeV, 1 TeV and 2 TeV, respectively. We can see from the figure that 
when the scale /i goes from 100 GeV to 1 TeV, the QCD X-factors vary from 1.176, 1.039, 0.900 to 
1.134, 1.029, 0.924 for ^ = 800 GeV, 1000 GeV, 2000 GeV, respectively. We can conclude that the 
theoretical uncertainty of the QCD correction at 0{as) level due to the variation of energy scale fj,, is 
under 4% for ^/s = 800, 1000, 2000 GeV, when rrih = 115 GeV and energy scale /x is in the range of 
100 GeV to 1 TeV. 

IV Summary 

In this paper we calculate the 0{as) QCD and ©(oew) one-loop electroweak radiative corrections to 
the e+e~ — > 77 — > tth^ process with Higgs boson in intermediate mass region at an e'^e~ linear collider 
(LC) in the SM. We investigate the dependence of the QCD and electroweak radiative corrections to 
both subprocess 77 tihP and parent process e^e~ ^ 77 ^ tihP on the Higgs boson mass mh 
and colliding energy \/l {^/s), and find that the QCD corrections can either increase or decrease the 
Born cross section, while the electroweak corrections always decrease the Born cross section of the 
e"'"e~ ^ 77 — > tth^ parent process and 77 — > tihP subprocess in the Higgs boson mass range 115 GeV< 
rrih < 200 GeV. We also notice that the 0{as) QCD corrections to process e~^e~ ^ 77 — > tth^ can be 
larger than the 0{aew) corrections depending on the Higgs boson mass rrih and e^e~ colliding energy 
y^. Both kinds of corrections may significantly decrease or increase the Born cross sections. The 
numerical results show that the 0{as) QCD relative corrections to the process e+e~ ^ 77 — > tth^ 
can reach 34.8% when ^/s = 800 GeV and mh = 200 GeV, while the O(aew) electroweak relative 



21 



corrections to the Born cross sections can reach —13.1%, —15.8% and —12.0% at ^/s = 800 GeV, 
1 TeV and 2 TeV, respectively. 

Acknowledgments: This work was supported in part by the National Natural Science Foundation 
of China and a grant from the University of Science and Technology of China. 

Appendix 

In Appendix we list the numerical comparison of the cross sections at the tree level for the process 
e+e^ ^ 77 ^ tth^ . In order to check our calculation we use two independently developed packages 
FeynArts 3|21j and ComDlIEPj22j to evaluate the cross sections. The results of ours and Kingman 
Cheung's are presented in Table 01 It is clear that our cross sections at tree level are not coincident 
with Cheung's results. 
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Figure Captions 

Figure^ The lowest order diagrams for the 77 tih^ subprocess.. 

Figure [2] The QCD pentagon diagrams for the 77 — > tth^ subprocess, whose ampHtudes include 
five-point tensor integrals of rank 4. The corresponding diagrams with interchange of the two incoming 
photons are not shown. 

Figure 121 The five-point pentagon electroweak one-loop diagrams for the 77 tihP subprocess, 
whose corresponding amplitudes include five-point tensor integrals of rank 4. 

Figure |3] The Born and one-loop QCD corrected cross sections for the 77 — > tihP subprocess as 
the functions of c.m.s. energy with m/j = 115, 150, 200 GeV, respectively. For each line type, 
the upper curve (in the energy region \/I > 1 TeV) is for the Born cross section and the lower one for 
the one-loop QCD corrected cross section. 

Figure [S] The Born and one-loop electroweak corrected cross sections for the 77 — > tth^ subprocess 
as the functions of c.m.s. energy (\/I) with m/j = 115, 150, 200 GeV. For each line type, the upper 
curve is for the Born cross section and the lower one for the one-loop electroweak corrected cross 
section. 
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Figure |51 The QCD one-loop relative corrections as the functions of c.m.s. energy for the 
77 tth^ subprocess with rrifi = 115, 150, 200 GeV, respectively. 

Figure ]7\ (a)The electroweak one-loop relative corrections as the functions of c.m.s. energy {^/I) 
with = 115, 150, 200 GeV for the 77 tih? subprocess. (b)The QED and weak one-loop relative 
corrections as the functions of c.m.s. energy with rrih = 115, 150, 200 GeV for the 77 — > tihP 
subprocess. For each line type, the upper curve is for the QED corrected cross section and the lower 
one for the weak corrected cross section. 

Figure El The QCD one-loop relative corrections as the functions of the mass of Higgs boson (m/^) 
for the 77 ^ tth^ subprocess with \/| = 500, 800, 1000, 2000 GeV, respectively. 

Figure iniThe electroweak one- loop relative corrections as the functions of the mass of Higgs boson 
(m/j) for the 77 tih? subprocess with ^/^ = 500, 800, 1000, 2000 GeV, respectively. 

Figure IIUI The Born and one-loop QCD corrected cross sections for the e"*"e~ — > 77 — > tthP 
process as the functions of c.m.s. energy {\fs) with nih = 115, 150, 200 GeV, respectively. For each 
line type, the upper curve(in the energy region ^/s > 1.4 TeV) is for the Born cross section and the 
lower one presents the one-loop corrected cross section. 

Figure Illl The Born and one-loop electroweak corrected cross sections for the e'^e~ ^ 77 ^ tthP 
process as the functions of c.m.s. energy {y/s) with rrih = 115, 150, 200 GeV, respectively. For each line 
type, the upper curve is for the Born cross section and the lower one presents the one-loop corrected 
cross section. 

Figure 1121 The QCD one-loop relative corrections to the e^e^ ^ 77 ^ tih^ process as the 
functions of c.m.s. energy {^/s) with nih = 115, 150, 200 GeV, respectively. 
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Figure ll3I The electroweak one- loop relative corrections to the e~^e~ — > 77 ^ tih^ process as the 
functions of c.m.s. energy {\/s) with mh = 115, 150, 200 GeV, respectively. 

Figure 1141 The QCD one-loop relative corrections to the e~^e~ — > 77 — > tth^ process as the 
functions of Higgs boson mass {rrih) with ^/s = 800, 1000, 2000 GeV. 

Figure 1151 The electroweak one-loop relative corrections as the functions of Higgs boson mass 
(mh) with ^ = 800, 1000, 2000 GeV, respectively. 

Figure ll6I The QCD relative corrections to the e~^e~ ^ 77 — > tth^ process as the function of the 
QCD renormahzation scale // , with ruh = 115 GeV and ^/s = 800, 1000, 2000 GeV. 
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Figure 3: The five-point pentagon electroweak one-loop diagrams for the 77 —>■ tth^ subprocess, whose 
corresponding amphtudes include five-point tensor integrals of rank 4. 
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Figure 4: The Born and one-loop QCD corrected cross sections for the 77 — > tih^ subprocess as the 
functions of c.m.s. energy {VS) with nih = 115, 150, 200 GeV, respectively. For each line type, the 
upper curve (in the energy region \/I > 1 TeV) is for the Born cross section and the lower one for the 
one-loop QCD corrected cross section. 



rrih (GeV) 


V§ (GeV) 


\S\max{%) 


115 


1000 


7.98 


130 


950 


9.15 


150 


890 


11.5 


170 


820 


14.0 


200 


710 


16.5 



Table 1: The maximum electroweak absolute relative corrections and the corresponding colliding 
energy positions for the 77 —>■ tthP subprocess with m/j = 115, 130, 150, 170, 200 GeV, respectively. 
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Figure 5: The Born and one-loop electroweak corrected cross sections for the 77 tthP subprocess as 
the functions of c.m.s. energy {\fs) with mh = 115, 150, 200 GeV. For each hne type, the upper curve 
is for the Born cross section and the lower one for the one-loop electroweak corrected cross section. 
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Figure 6: The QCD one-loop relative corrections as the functions of c.m.s. energy {yfs) for the 
77 — tthP subprocess with ruh = 115, 150, 200 GeV, respectively. 
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Figure 7: (a)The electroweak one-loop relative corrections as the functions of c.m.s. energy with 
mh = 115, 150, 200 GeV for the 77 — tth^ subprocess. (b)The QED and weak one-loop relative 
corrections as the functions of c.m.s. energy (\/^) with rrifi = 115, 150, 200 GcV for the 77 tth^ 
subprocess. For each line type, the upper curve is for the QED corrected cross section and the lower 
one for the weak corrected cross section. 



32 



1 1 


1 1 

^Ps = 500 GeV 




^fi = 800 GeV _ 




/V -I C\C\C\ T 7 

vs — 1000 GeV 




\/l = 2000 GeV _ 


" 1 r ' 


1 1 



100 120 140 160 180 200 

mh [GeV] 

Figure 8: The QCD relative corrections to tire cross sections of 77 tthP subprocess as the functions 
of the mass of Higgs boson {rrih) with y/§ = 500, 800, 1000, 2000 GeV, respectively. 
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Figure 9: The electroweak relative corrections to the cross sections of 77 tthP subprocess as the 
functions of the mass of Higgs boson (m^) with V§ = 500, 800, 1000, 2000 GeV, respectively. 
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Figure 10: The Born and one-loop QCD corrected cross sections for the e+e" ^ tth^ process 

as the functions of c.m.s. energy (y^) with nih = 115, 150, 200 GeV, respectively. For each line type, 
the upper curve(in the energy region ^/s > 1.4 TeV) is for the Born cross section and the lower one 
presents the one-loop corrected cross section. 



mi, (GcM) 


^ (Gi-V) 


\(^\nu,r(%) 


115 


2000 


7.10 


130 


1700 


8.28 


150 


1435 


10.7 


170 


1220 


13.4 


200 


1020 


15.9 



Table 2: The maximum absolute relative corrections and the corresponding colliding energy y/s posi- 
tions for the e+e~ — 77 — tth*^ process with m/j = 115, 130, 150, 170, 200 GeV, respectively. 
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Figure 11: The Born and one- loop electroweak corrected cross sections for the e^e" ^ 77 ti}^ 
process as the functions of c.m.s. energy {\fs) with m/i = 115, 150, 200 GeV, respectively. For each line 
type, the upper curve is for the Born cross section and the lower one presents the one-loop corrected 
cross section. 
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Figure 12: The QCD relative corrections to the e^e ^ 77 ^ tt}{^ process as the functions of c.m.s. 
energy {\fs) with mh = 115, 150, 200 GeV, respectively. 
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Figure 13: The electroweak relative corrections to the e^e — ^ 77 — ^ ttY^ process as the functions of 
c.m.s. energy {\fs) with = 115, 150, 200 GeV, respectively. 
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Figure 14: The QCD relative corrections to the e+e — >^ 77 — >^ ttY^ process as the functions of 
Higgs-boson mass (m^) with -y/s = 800, 1000, 2000 GeV, respectively. 
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Figure 15: The electrowcak relative corrections to the e+e ^ 77 ^ tth^ process as the functions of 
Higgs-boson mass (mh) with -s/s = 800, 1000, 2000 GeV, respectively. 
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Figure 16: The QCD relative corrections to the e'^e~ ^ tthP process as the function of the 

QCD renormalization scale /x with = 115 GeV and y/s = 800, 1000, 2000 GeV. 
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mt [GeV] 


mh [GeV] 


yfs [GeV] 


(7 [fb] (Ref.jin]) 


a [fb] (FeynArts) 


a [fb] (CompHEP) 


120 


60 


500 


0.45 


0.391(0) 


0.391(8) 






1000 


2.6 


2.18(7) 


2.19(1) 






2000 


2.8 


2.39(1) 


2.39(1) 


150 


60 


1000 


3.2 


2.74(1) 


2.74(5) 






2000 


4.1 


3.42(1) 


3.42(2) 




140 


1000 


0.36 


0.311(8) 


0.311(6) 






2000 


0.95 


0.805(9) 


0.805(6) 


180 


140 


1000 


0.40 


0.341(3) 


0.341(5) 






2000 


1.2 


1.05(5) 


1.05(5) 



Table 3: The numerical comparison of the cross sections of the process e^e ^ 77 ^ tth^ at tree-level 
with the results in Ref . |20j by using FeynArts 3 and CompHEP packages. 
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